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ABSTRACT 

Here we present the colours of individual globular clusters (GCs) around the 
massive elliptical galaxy NGC 1407 out to a projected galactocentric radius of 140 
kpc or 17 galaxy effective radii (Re)- Such data are a proxy for the halo metallicity. 
We find steep, and similar, metallicity gradients of ~ -0.4 dex per dex for both the blue 
(metal-poor) and red (metal-rich) GC subpopulations within 5-8.5 Re (40-70 kpc). At 
larger radii the mean GC colours (metallicity) are constant. A similar behaviour is seen 
in a wide-field study of M87's GC system, and in our own Galaxy. We interpret these 
radial metallicity trends to indicate an inner region formed by early in-situ dissipative 
processes and an outer halo formed by ongoing accretion of low mass galaxies and their 
GCs. These results provide observational support for the model of galaxy formation 
whereby massive galaxies form inside-out in two phases. We have also searched the 
literature for other massive early-type galaxies with reported GC metallicity gradients 
in their inner regions. No obvious correlation with galaxy mass or environment is found 
but the sample is currently small. 

Key v^rords: globular clusters: general - galaxies: star clusters - galaxies: individual 
(NGC 1407) - galaxies: formation 



1 INTRODUCTION 

An emerging picture of galaxy formation is that mas- 
sive early-type galaxies grow from relatively compact 'seed' 
galaxies with effective radii of ~ 1 kpc at z ^ 2 (sometimes 
referred to as 'red nuggets') to today's massive systems of 
effective radii of ~ 5 kpc (Oser et al. 2010; Gobat et al. 
2010). This size evolution has a redshift dependence of the 
form (1-l-z)"^-^ (Franx et al. 2008; Williams et al. 2010) and 
is driven by the growth of the outer regions of galaxies that 
experience ongoing accretion (van Dokkum et al. 2010). 

According to theory, the growth of stellar mass occurs 
in two phases: 1) An m-situ dissipative phase with cold flows 
(Keres et al. 2005; Dekel et al. 2009; Naab et al. 2009; Oser 
et al. 2010) that is reminiscent of a 'monolithic collapse' 
(Larsen 1975) and 2) ongoing accretion from minor merg- 
ers of lower mass systems (Abadi et al. 2006; Hopkins et al. 
2009; Zolotov et al. 2010). The in-situ stars form, and are 
in place, early in the Universe (i.e. ^ 10 Gyrs ago), whereas 
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the outer halo experiences a more continous growth as addi- 
tional material is accreted over time. The transition between 
the m-situ dominated inner regions and the accretion dom- 
inated outer halo is predicted to occur at radii of tens of 
kiloparsecs (Abadi et al. 2006; Kolotov et al. 2010), with 
more massive galaxies undergoing more accretion and hav- 
ing larger transition radii (Oser et al. 2010). 

This theoretical picture is difficult to test observation- 
ally. For very nearby (~ 10 Mpc) early-type galaxies, indi- 
vidual stars can be resolved by the Hubble Space Telescope 
(HST) and the radial metallicity distribution probed. For 
example, in NGC 3379 the halo shows a metallicity transi- 
tion at 10-13 Re, becoming metal-poor (Harris et al. 2007). 
Only a handful of massive galaxies can be studied in this 
way. 

Intergrated light observations are required to obtain 
the metallicity distribution for more distant galaxies. How- 
ever, the surface brightness profiles of early-type galaxies fall 
off rapidly with galactocentric radius making it difficult to 
probe beyond a few effective radii. An example is the recent 
work of Coccato et al. (2010) who, after 13.5 hrs integration 
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on an 8m telescope, obtained stellar population gradients 
in the massive Coma elliptical galaxy NGC 4889 to an un- 
precedented 4 Re. They derived metallicity, alpha element 
and age gradients from an analysis of the absorption lines 
using Lick indices. They found a "break" radius at ^1.2Ke 
(18 kpc) dividing the inner galaxy that exhibited a steep 
metallicity gradient and high alpha element ratio (indica^ 
tive of rapid star formation) from the outer region that had 
a metallicity gradient consistent with zero and evidence for 
more prolonged star formation. They described the inner re- 
gion as formed in an early dissipative collapse and the outer 
region as having an accreted origin. Given the telescope time 
involved it will be difficult to repeat such a study for a large 
sample of ellipticals. 

Another tracer of halo stellar populations is therefore 
required. Globular clusters (GCs) have many advantages in 
this regard. As some GCs are destroyed over time, they con- 
tribute directly to the halo field star population. In the case 
of the Milky Way, about 1/3 of the original GC system may 
have been destroyed (Mackey & van den Bergh 2005) mak- 
ing up half of the current halo stellar mass (Martell & Grcbel 
2010). When a low mass galaxy is accreted by a larger one, 
GCs can survive the accretion process (Forbes & Bridges 
2010) and the properties of these GCs contain information 
about their original host galaxy. Also, GC systems can ex- 
tend beyond 100 kpc and hence they probe galaxy halos to 
many effective radii. 

The globular cluster systems of large galaxies are gener- 
ally found to consist of two subpopulations - blue (or metal- 
poor) and red (or metal-rich) . These are thought to be asso- 
ciated with the halo and bulge/spheroid component of galax- 
ies, respectively (Brodie & Stradcr 2006). The mean colour 
of both subpopulations have been shown to correlate with 
the luminosity or mass of the host galaxy, with some possi- 
ble flattening of the relation for low mass galaxies (Forbes, 
Brodie & Grillmair 1997; Larsen et al. 2001; Strader, Brodie 
& Forbes 2004; Peng et al. 2006). We note that the mean 
colour used in such studies is generally derived from the 
central regions of a galaxy. Recent studies of extragalactic 
GC systems have revealed two additional trends: the colour 
(metallicity) of GCs with their magnitude (mass) and colour 
with projected galactocentric radius. 

The first trend, called the 'blue tilt' (a tendency for the 
blue GCs to have redder colours at brighter magnitudes), 
was first discovered using the Advanced Camera for Surveys 
(ACS) on the HST to examine extragalactic GC systems. 
The ACS provides highly accurate photometry and GC can- 
didate lists that are largely free of contaminants (due to the 
fact that GCs are partially resolved out to distances of ~ 
20 Mpc). The blue tilt has now been reported in a vari- 
ety of galaxies including the most massive ellipticals (Harris 
et al. 2006; Harris 2009a), lower mass ellipticals (Strader 
et al. 2006), early- type dwarfs (Mieske et al. 2006a), early- 
type spirals (Spitler et al. 2006) and a Milky Way like spiral 
(Forbes et al. 2010). It has not been detected in the Milky 
Way GC system (most likely due to the small number of 
massive GCs) nor the massive elliptical NGC 4472 (Strader 
et al. 2006; Mieske et al. 2006a) despite its large GC system. 

Given the evidence that extragalactic GCs are mostly 
old (Brodie & Strader 2006), the blue tilt implies a mass- 
metallicity relation for the metal-poor subpopulation of 
GCs. Various explanations have been proposed (see Bekki 



et al. 2007; Mieske et al. 2006a; Mieske 2009), but perhaps 
the most plausible explanation to date is one of self enrich- 
ment, whereby more massive GCs are enriched with heavier 
metals during their brief formation period (Parmentier 2004; 
Strader & Smith 2008; Bailin & Harris 2009). 

The second trend, whereby the mean colour of each GC 
subpopulation varies with galactocentric radius, was first 
seen in ground-based photometry of M49 using metallicity- 
sensitive filters (Geisler et al. 1996). Such a trend needs to be 
carefully separated from the trend of a declining mean colour 
with radius in the overall GC system (due to the changing 
relative proportions of blue and red subpopulations), and 
has onlj' been measured for a handful of other systems (Forte 
et al. 2001; Bassino et al. 2006; Harris 2009a,b). These stud- 
ies indicate a negative metallicity gradient for both the blue 
and red subpopulations in elliptical galaxies. In the case of 
the Milky Way, a metallicity gradient is seen in both GC sub- 
populations within the central 10 kpc of slope « -0.30 dex 
per dex (Harris 2001). Beyond 10 kpc, the mean GC metal- 
licity in the Galactic halo is constant at [Fc/H] ~ -1.5. The 
M31 GC system reveals a strong metallicity gradient within 
50 kpc, but flattens out to a near constant metallicity of 
[Fe/H] ~ 1.6 (Alves-Brito et al. 2009). These strong metal- 
licity gradients in the inner halo regions suggest a dissipative 
formation process for both GC subpopulations. 

Here we present the radial distribution of GC colour 
(and hence metallicity) using deep wide-field Subaru imag- 
ing of the massive elliptical galaxy NGC 1407 out to ~140 
kpc or 17 Re. We find a strong radial gradient which is sup- 
ported by HST imaging and Keck spectroscopy. We com- 
pare the radial gradient to those of other massive ellipticals, 
finding that the gradient in the NGC 1407 GC system is 
the steepest published to date. We briefly discuss the wider 
implications for galaxy formation from our results. 



2 NGC 1407 

NGC 1407 is the central massive elliptical in a large group 
(Brough et al. 2006) with a rich globular cluster system 
(Forbes et al. 2006). The galaxy itself reveals a steep metal- 
licity gradient in the starlight (Spolaor et al. 2008; Foster et 
al. 2009) and is uniformly old with a supersolar alpha ele- 
ment ratio within one effective radius. Our previous optical 
spectra of 20 bright GCs in NGC 1407 indicated that they 
are mostly old (~11 Gyr) with a range of metallicity -1.5 
< [Z/H] < 0.0 (Cenarro et al. 2007). Some of the spectra 
showed evidence for the presence of blue horizontal branch 
stars. Our previous dynamical study of the GC system (Ro- 
manowsky et al. 2009) found weak evidence for rotation and 
a bias towards tangential orbits in the outer regions. We de- 
rived a total virial halo mass of 6 x 10^^ Mq which indicates 
that NGC 1407 and its surrounding group is extremely dark 
(with M/Ls ~ 800). We assume an effective radius for NGC 
1407 itself of 1.17 arcmin as adopted by Harris (2009) and 
Foster et al. (2009). 

Images of NGC 1407 have been taken using the ACS 
camera on board the HST using F435W (B) and F814W (I) 
filters covering 3.5 arcmin x 3.5 arcmin. An analysis of the 
GC system using these data has been reported by Forbes et 
al. (2006) and Harris et al. (2006) and Harris (2009a). The 
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GC system reveals a classic blue/red bimodality with a blue 
tilt in the blue subpopulation. 

Harris (2009a) studied the colour-radius distribution for 
the NGC 1407 GC system finding an insignificant (la) trend 
for both the blue and red GCs examined. However, when 
the NGC 1407 GC system was combined with that of 5 
other massive elliptical galaxies, a statistically significant 
metallicity gradient of -0.10 dex per dex was seen at >4.5cr 
for both subpopulations. 

The distance to NGC 1407 is somewhat uncertain. Here 
we adopt the mean value of 15 distances given by NED, i.e. 
m-M = 31.85 ± 0.37 (or 23.77 ± 3.8 Mpc). At this distance, 
1 arcmin equals 6.9 kpc. 



3 PHOTOMETRIC DATA 

Here wc use wide-field (34 arcmin x 27 arcmin) g,r,i im- 
ages taken using the Suprime-Cam camera on the Subaru 
8m telescope. Total exposure times were 3hr, Ihr, 0.9hr for 
the g,r,i images respectively. The seeing was good at around 
0.5 arcsec in each filter. The Suprime-Cam has a pixel scale 
of 0.2 per pixel. Full details of the observations, data re- 
duction and analysis will be presented in Spitler et al. (2011, 
in prep.). Here we give a brief summary. 

The images were reduced and combined using the 
Suprime-cam standard pipeline software. The final mosaic 
was placed on the USNO-B2 astrometric system using com- 
mon stars. Before the initial object detection, a smooth 
model of NGC 1407 (and the galaxy NGC 1400, which lies at 
a projected separation of 12 arcmin) was subtracted, using 
the ELLIPSE task in IRAF to help with object detection. 
DAOFIND was then used to select objects with a threshold 
^3 times the standard deviation of the background count 
level. 

For each point-like object an aperture magnitude and 
position was measured. Aperture magnitudes were corrected 
to total magnitudes assuming that the objects were point 
sources (a reasonable assumption for GCs with an average 
expected size of 3 pc or 0.03 arcsec under 0.5 arcsec seeing) 
and applying filter-dependent aperture corrections. These 
were converted into standard magnitudes using zeropoints 
determined with standard star observations. Extinction cor- 
rections using the maps of Schlegel et al. (1998) were ap- 
plied assuming no variation across the image as NGC 1407 
is known to be relatively free of dust from mid-lR obser- 
vations (Temi et al. 2009). Hereafter we quote extinction 
corrected magnitudes and colours. A colour selection was 
applied so that objects within la of the GC sequence in g-r 
and r-i colour-colour parameter space (i.e. 0.4 < g-i < 1.4 
and 0.3 < g r < 0.9) were included. 

Finally, any GC candidates associated within 6 arcmin 
radius of NGC 1400 were removed from the object list for the 
subsequent analysis. This effectively removes any possible 
contamination from the GC system of NGC 1400. 



4 RESULTS 

In Figure 1 we show the distribution of GC colour with 
projected galactocentric radius out to 20 arcmin from our 
wide-field Subaru imaging. As well as the individual data 



points we show the mean value and error on the mean in 
several radial bins using the NMIX statistical test (Richard- 
son & Green 1997). The figure shows only those GCs with 
i < 24 (Mi > -7.9). The data show evidence for a clear 
negative metallicity gradient out to ~6-10 arcmin (5-8. 5Re 
or 40-70 kpc) and are consistent with zero gradient beyond 
that. A simple linear fit to radii within 10 axcmin gives cen- 
tral colours of g-i = 0.880 ± 0.011 and 1.186 ± 0.007 for 
the blue and red subpopulations respectively. Similarly, the 
slopes arc found to be 0.0084 ± 0.022 and 0.013 ± 0.002 
mag per arcmin. Thus the blue subpopulation gradient is at 
the ~4iT level and the red one at the ~6cr. 

Following the method of Harris (2009b) we also fit the 
entire colour gradient with the form g-i = a + b log(R/Re), 
where Re — 1.17 arcmin. We fit only the data points interior 
to ~10 arcmin (although including radial bins beyond this 
radius makes little difference to the final values). For the 
blue subpopulation we find a = 0.818 ± 0.01 and b — -0.063 
± 0.01 mag per dex. Similarly for the red subpopulation we 
find a = 1.12 ± 0.02 and b = -0.068 ± 0.02 mag per dex. For 
both subpopulations the uncertainty quoted comes from the 
range of fits that include/exclude the innermost data point 
(the steeper slope being associated with an exclusion of the 
innermost radial bin). Within the combined errors, the blue 
and red GC subpopulations show a similar colour gradient. 

To convert the g-i colour gradients into metallicity gra- 
dients we follow Harris (2009a), i.e. Ag-i = 0.21A[Fe/H], 
giving metallicity gradients of -0.38 ± 0.06 dex per dex for 
the blues and -0.43 ± 0.07 dex per dex for the reds. 

We also show in Figure 1 metallicities corresponding to 
g-i colours assuming the linear relation found by Sinnott et 
al. (2010). This relation is in part based on Milky Way GC 
metallicities from Harris (1996) which is close to the Car- 
retta & Gratton (1997) scale and hence we assume is closer 
to total metallicity [Z/H] (see also discussion by Mendel et 
al. 2007). 

We have calculated the probability density function 
(PDF) for the data shown in Figure 1 in six radial bins (cen- 
tred on 0.92, 1.81, 2.72, 3.85, 5.23 and 7.04 arcmin) of equal 
GC number. The PDF is sensitive to colour substructure 
and is shown graphical^ in Figure 2. Other than the inner- 
most radial bin, the PDF for the blue GCs confirms their 
radial colour gradient. For the red GCs the radial trend is 
less clear due to possible colour substructure in the inner 
regions. 

In Figure 3 wo show the distribution of GC colour with 
galactocentric radius from the HST ACS imaging (Forbes 
et al. 2006). These data have the advantage of having a 
lower contamination rate than the Subaru imaging (which 
we estimate to be less than 5% for bright GCs based on 
follow-up spectra; Romanowsky et al. 2009) but probe to 
much smaller galactocentric radii. Again we show the mean 
colours and the error on the mean in various radial bins 
with equal numbers of GCs. A clear colour gradient is seen 
in the red subpopulation and to a lesser extent in the blue 
subpopulation supporting our Subaru results above. Fitting 
B-1 = o -I- 6 log(R/Re) we find a = 1.62 ± 0.02 and b = 
-0.147 ± 0.065 dex per dex for the blue subpopulation and 
a = 2.05 ± 0.01 and 6 = -0.169 ± 0.005 dex per dex for the 
red GCs. 

Assuming a conversion of B-I colour into metallicity 
of AB-I = 0.375A[Fe/H] (Harris 2009b), these slopes give a 
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mctallicity gradient of -0.39 ± 0.17 dex per dex for the blues 
and -0.45 ± 0.01 for the reds. These values are identical to 
those inferred from the larger radii g-i gradient. 

We also show in Figure 3 the metallicities correspond- 
ing to B-I colours assuming the linear relation quoted by 
Harris (2009b). As for the g-i colours, this relation based 
on Milky Way GC metallicities from Harris (1996) which is 
close to the Carretta & Gratton (1997) scale and hence total 
metallicity [Z/H]. 

The radial stellar population properties for the NGC 
1407 field stars out to ~0.5Re have been measured, from 
an analysis of the absorption lines, by Spolaor et al. (2008). 
They found the galaxy to have a central total mctallicity of 
[Z/H] ~ +0.3 dex and a steep mctallicity gradient of -0.38 
± 0.04 dex per dex. Those values are similiar to those found 
for the red GCs in NGC 1407 and are illustrated in Figure 
3. The error on individual points from Spolaor et al. is ± 
0.15 dex. This supports previous claims that the red GCs 
are closely associated with the bulge/spheroid component 
of elliptical galaxies (Forbes & Forte 2001; Forte et al. 2009; 
Spitler 2010). The metalhcity gradient of Foster et al. (2009) 
has a similar slope but has individual errors of ± 1 dex (and 
is not shown). 

Could the colour gradient for the blue subpopulation 
seen in Figures 1, 2 and 3 simply be due to a changing 
strength of the blue tilt with radius? For this to be the 
case the blue tilt would need to be significantly stronger 
(i.e. tilted more to the red at bright magnitudes) in the 
inner galaxy regions and to disappear (i.e. no blue tilt) for 
galactocentric radii beyond ~10 arcmin. Mieske et al. (2006, 
2010) found weak evidence for a stronger blue tilt at small 
galactocentric radii in their data for Virgo and Fornax early- 
type galaxies. However, for the large galaxies in their study, 
they were generally restricted to the inner ~1.5Re due to the 
HST ACS field-of-view. Harris (2009a) explored this issue 
out to ~5Re using a composite of the GC systems in 6 giant 
ellipticals, finding no convincing evidence for a stronger blue 
tilt in the inner regions. Further evidence comes from the 
fact that NGC 4472 reveals a radial gradient (Geisler et al. 
1996) but no blue tilt (Strader et al. 2006). 

We have divided the Subaru data into 3 radial bins, i.e. 
to a projected galactocentric radius of 5 arcmin, between 5 
arcmin and 10 arcmin and beyond 10 arcmin. All three ra- 
dial bins reveal a blue tilt of similar strength. The relatively 
constant strength of the blue tilt with radius can not explain 
the observed gradient in the mean colour with radius within 
10 arcmin (nor the lack of a gradient beyond 10 arcmin). We 
thus conclude that the colour gradient observed in the blue 
subpopulation is not caused by a changing colour-magnitude 
trend with radius. Finally, the radial metallicity gradient of 
the blue subpopulation is similar to that of the red GCs, 
which do not show a 'red tilt'. 

Could the GC colour gradients be due to contamina- 
tion? The radial surface density profiles of the blue and red 
GC candidates around NGC 1407 (Spitler et al. 2011, in 
prep), show that they are dominated by objects that are as- 
sociated with NGC 1407 for radii <10 arcmin (i.e. for which 
the strong metallicity gradient is seen). Based on these pro- 
files the estimated contamination rate at 10 arcmin is less 
than 20% for the blue GCs but could be as high as 50% for 
the red GCs. Thus the transition radius for the blue GCs 
appears to be robust but the flattening of the gradient for 



the red GCs at large radii could be an artifact of increased 
contamination. 

Foster et al. (2010) explored the use of the Calcium 
Triplet (CaT) lines to derive spectroscopic metallicities for 
GCs in NGC 1407 but did not explore the radial variation of 
CaT metallicity for the GC subpopulations. Using a linear 
CaT to metallicity transformation based on Milky Way GCs 
by Armandroff & Zinn (1988), they found the red GCs to 
have reasonable metallicities but the blue GCs had inferred 
metallicities that were too high, i.e. they had stronger CaT 
strengths than expected. This could be due to an extra con- 
tribution from blue horizontal branch stars (as suggested 
by Cenarro et al. 2007 for the NGC 1407 GCs). In Figure 
4 we plot the CaT derived metallicities (converted into to- 
tal metallicity assuming [Z/H] — [Fe/H] -I- [a/Fe], where 
[a/Fe] is assumed to have a constant value of 0.3 typical of 
Milky Way GCs). The data, taken with the DEIMOS spec- 
trograph on the Keck telescope, have been restricted to GCs 
with S/N ^12. We show the red and blue GCs separately. 
The blue GCs show significant scatter about a mean value of 
[Z/H] ~ -0.5, again suggesting a problem with their derived 
metallicities (see Foster et al. 2010 for a full discussion for 
this and possible explanations). The red GCs, on the other 
hand, show a clear radial gradient. Here we overplot the total 
metallicity gradient derived from Figure 1 on the red GCs. 
Although not a fit to the data, the metallicity gradient based 
on g-i colour provides a remarkably good representation of 
the overall red GC trend. 



5 DISCUSSION 

From deep wide-field imaging we find a transition from a 
strong colour gradient in both GC subpopulations which 
flattens to a constant value, with the transition occurring at 
a projected radius of around 6-10 arcmin. Approximately 
half of the GCs are located within this radius. At a slightly 
larger radius (~12 arcmin) Zhang et al. (2007) report that 
the X-ray surface brightness around NGC 1407 declines to 
a near constant value, which they refer to as the 'gas halo 
boundary'. 

The behaviour of the NGC 1407 GC colour profiles are 
similar to those in the Milky Way, for which strong gradi- 
ents (-0.30 dex per dex; Harris 2001) are seen in both the red 
and blue GCs interior to 10 kpc and a constant metallicity 
of [Fe/H] ~ -1.5 for the blue GCs at larger radii (there are 
essentially no red GCs beyond 10 kpc). The gradient seen 
in the NGC 1407 GC system is the strongest reported to 
date in an elliptical galaxy, with the previous strongest gra- 
dient being -0.28 dex per dex in the NGC 1399 GC system 
reported by Bassino et al. (2006). We note that actual 3D 
gradients may be stronger than those observed in projection. 

Most previous detections of GC colour gradients in 
early-type galaxies have been generally restricted to ^ 5Re 
and may not extend far enough to probe the accretion dom- 
inated region of the halos of massive galaxies. One excep- 
tion is the wide-field study of M87 by Harris (2009a) using 
CFHT Megacam which allowed the GC system to be probed 
to ~20Re. To quote Harris This [metallicity gradient] feature 
is clearest for R 8 Re or about 60 kpc; beyond that point, 
the red clusters are almost absent and the blue clusters show 
no significant change in mean color. Thus the GC system of 
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M87 also appears to reveal a transition, at a similar scale in 
both eflfective radius and kpc to the NGC 1407 GC system. 

The transition from a strong GC metallicity gradient 
to a constant value may represent a transition from a dis- 
sipativc process (e.g. a gaseous collapse or gas-fed flows) to 
a dissipationless one (e.g. accretion of low mass galaxies). 
This is in qualitative agreement with the recent simulations 
of galajcy formation in two phases (e.g. Naab et al. 2009; 
Oser et al. 2010). In this model, the first in-situ phase is as- 
sociated with a dissipative collapse, which is followed by an 
extended phase of ongoing accretion of low mass galaxies. 
More massive galaxies have larger transition radii. In this 
picture galaxies are built from the inside-out. In the dissipa- 
tional merger simulations of Hopkins et al. (2009) the outer 
constant metallicity is determined by the metallicity of the 
progenitor galaxies. 

An alternative explanation is that the observed colour 
profile is caused by a changing relative mix of inner GCs, 
that have fixed colours, and outer GCs that are bluer on av- 
erage. Thus as the fraction of outer bluer GCs increases with 
radius a mean colour gradient is observed until the transition 
radius is reached, when the outer GCs completely dominate 
the overall GC colours (beyond 70 kpc in the case of NGC 
1407). Such a situation could in principle be setup by a se- 
ries of dry mergers without invoking dissipation, however 
it would need to be fine-tuned in the sense of giving simi- 
lar observed metallicity gradients for both the blue and red 
subpopulations and to occur at similar transition radii. 

A kinematic analysis of 172 GCs around NGC 1407 
revealed a bias towards tangential orbits with increasing 
galactocentric radius (Romanowsky et al. 2009). This or- 
bital anisotropy is difficult to reconcile with either mergers 
or an accretion origin for the GCs. However, the study was 
largely confined to radii within 8 arcmin. 

We note that kinematic transitions at large radii have 
now been reported in handful of massive early-type galaxies 
based on analysis of their stellar light and planetary nebulae 
as well as their GC systems (Proctor ct al. 2009; Coccato 
et al. 2009; Foster et al. 2010; Arnold et al. 2010). However, 
more work is needed to increase radial sampling and number 
of galaxies studied. 

Prompted by trends between the metallicity gradient 
of the galaxy starlight and galaxy mass from both observa- 
tions (e.g. Spolaor et al. 2009) and theory (e.g. Tortora et 
al. 2010), we have collated the statistically significant (i.e. 
more than 2a) GC metallicity gradients from the literature 
and central velocity dispersions (a proxy for galaxy mass). 
Table 1 shows that there is no obvious trend between the 
GC metallicity gradient and the mass of the host galaxy (as 
traced by the central velocity dispersion), however the sam- 
ple is small and covers a limited range in host galaxy mass 
and environment. 

The presence of two GC subpopulations in most galax- 
ies is strong evidence that galaxies did not form by a simple 
monolithic collapse. The two phases of formation for NGC 
1407 described here, support this view. Indeed, the in-situ 
GCs appear to have formed in an early dissipative event as 
envisaged by Forbes et al. (1997). Here we have shown that 
both GC subpopulations appear to have taken part giving 
rise to similar strength metallicity gradients and transition 
radii. The accreted GCs provide an important contribution 
as envisaged in the scenario of Cote et al. (1998). However in 



their scenario for the GC systems of massive galaxies, only 
the red GCs were involved in a dissipative formation event 
and the blue GCs were all accreted. 



6 CONCLUSIONS 

Using deep wide-field imaging from the Subaru Suprime- 
cam we have explored the variation of GC colour (and hence 
metallicity) with galactocentric radius. We detect GCs out 
to -^20 arcmin or 17 galaxy effective radii (Re)- This corre- 
sponds to ~ 140 kpc for an assumed distance of 23.77 Mpc. 
This is therefore one of the most extended radial studies of 
a GC system to date. 

We find that both the blue and red GC subpopulations 
reveal a strong colour gradient within ~6-10 arcmin (40-70 
kpc) or 5-8.5 Re. Beyond this, the mean GC colours are con- 
stant. The presence of a colour gradient within the inner 3 
arcmin is supported by HST ACS imaging. When converted 
into metallicity, the gradient for the red GC subpopulation 
is similar to that for the field stars in the central ~0.5Re. We 
also report a spectroscopic metallicity gradient in the inner 
10 arcmin based on Calcium Triplet measurements from the 
Keck DEIMOS spectrograph. 

The change from a strong GC metallicity gradient to 
a constant metallicity may indicate the transition between 
an inner halo formed by in-situ dissipative processes, with 
both GC subpopulations taking part, to an outer one that 
is largely due to accretion of minor mergers. This transi- 
tion radius is similar to the hot gas boundary seen in X-ray 
brightness profile of NGC 1407. 

We find that the colour gradient seen in the blue GC 
subpopulation is not due to a changing 'blue tilt' strength 
with radius. Therefore, if the blue tilt is due to a self- 
enrichment process it does not appear to have a strong ra- 
dial, or environmental, dependence. 

We have searched the literature for statistically signif- 
icant GC metallicity gradients in elliptical galaxies. We do 
not find a trend with host galaxy mass; however the sample 
is small and largely restricted to massive galaxies in clus- 
ters/groups. We note that previous work on the correla- 
tion between mean GC colour and host galaxy mass has 
generally been restricted to galaxy central regions. Future 
wide-field studies of GC system radial colour gradients in a 
wider range of host galaxies and environments may reveal 
new trends and hence further insight into galaxy formation 
process. Such a study is relatively inexpensive in terms of 
telescope time. 

Our results provide observational support for the mass 
growth of massive elliptical galaxies in two phases, i.e. early 
dissipative collapse followed by ongoing accretion of low 
mass systems (e.g. Oser et al. 2010). In other words, el- 
liptical galaxies form inside-out. 
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Table 1. Globular cluster system metallicity-radius gradients in units of dex per dex. The central velocity dispersion (in km s~^) comes 
from the Hypcrlcda database. The source of each globular metallicity gradient is given in the last column. The metallicity gradients are 
converted from colour gradients assuming an old stellar population. The la uncertainty when available is given in brackets. 



Galaxy 


Environ. 




Blue 


Red 


Source 


NGC 1407 


Group 


272 


-0.38 (0.06) 
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Geisler et al. (1996) 
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Forte et al. (2001) 
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Harris (2009a) 
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Figure 1. Globular cluster colour-radial distribution from Subaru photometry. The small circles show individual GCs from Subaru 
photometry with i < 24 (and errors ^ 0.03 mag). Mean colours and error on the mean are represented by filled squares. The solid lines 
show the best fit linear gradient within 10 arcmin. The top and right axes show transformations into effective radius and metallicity 
respectively. Both subpopulations show a strong gradient in mean colour within the central ~10 arcmin (8.5 Re, 69 kpc). The galaxy 
starlight at 1 Re has a colour g— i ~ 1.2, similar to the mean colour of the red GCs. 
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Figure 4. Globular cluster metallicity-radial distribution derived from calcium triplet measurements. The top panel shows the red GC 
subpopulation, and the bottom panel shows the blue GC subpopulation. A clear radial gradient is seen in the red GC subpopulation. 
The solid line shown in the top panel is not a fit to the data but rather the metalicity gradient from Figure 1 superposed. 



